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Abstract 

The CLEO Collaboration has made the first observations of hadronic transitions among bot- 
tomonium ihh) states other than the dipion transitions among T(nS) states. In our study of 
T(3S) decays, we find a significant signal for T(3S)^7a;T(lS) that is consistent with radiative 
decays T(3S)^7X6i.2(2i-'), followed by X6i,2(2P)— >a;T(lS). The branching ratios we obtain are 
S(xw(2P)-^T(lS)) = (1.63t[J:itO:}i)% and ^?(x62(2P)^a;T(lS)) = (1.10l°;ii°;li)%, in which 
the first error is statistical and the second is systematic. 
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The only hadronic decays of bottomonia that have been experimentally observed to date 
are the tttt (tttt = it~^it~ and 7r°7r°) transitions among the T(nS) states [if. In Fig. [Hwe 
show the spectrum of bottomonium states below open bottom threshold. 

Hadronic transitions among heavy quarkonia are generally understood to proceed by the 
emission of soft gluons, and subsequent hadronization of the gluons. The analysis of heavy 
quarkonium hadronic transitions is one of few possible laboratories for the study of the 
physics of the soft gluon emission and hadronization process that governs such decays. 

Most theoretical work dedicated to these transitions has been built around a multipole 
expansion of the color field, an idea first proposed by Gottfried and Yan 0. A fairly 
substantial literature exists which attempts to describe in detail the tttt transitions that 
have been observed While these transitions do provide important information about 
strong interaction dynamics in heavy quark systems, the investigation of other hadronic 
decay modes {i.e., involving r], cu or multiple vr) should offer a different perspective. 

In the multipole expansion model, a hadronic transition involving uj requires three gluons 
in an El * El * El configuration For such a purely electric coupling, Voloshin jsj has 
recently predicted roughly equal rates for the decay of the two states Xbii'^P) and Xb2(2-P) to 
a;T(lS). In this Letter, we report on the observation of the transitions xw,2(2P)— >co'T(lS). 

The data set consists of (5.81 ± 0.12) x 106T(3S) decays observed with the CLEO III 
detector at the Cornell Electron Storage Ring. Charged particle tracking is done by the 
47-layer drift chamber and a four-layer silicon tracker which reside in a 1.5 T solenoidal 
magnetic field. Photons are detected using an electromagnetic calorimeter consisting of 
7784 CsI(Tl) crystals distributed in a projective barrel geometry. The particle-identification 
capabilities of the CLEO III detector are not used in the present analysis. 

Events consistent with the final state of 77r"''7r~7r°£'*"£~ are selected by requiring two 
high momentum charged tracks (p > 4 GeV/c) and two or three charged low momentum 
tracks (0.12 < p < 0.75 GeV/c). The low momentum tracks are required to come from the 
interaction region using criteria obtained by studying charged pion tracks in a sample of 
events from the kinematically similar decay T(2S)^7r+7r~T(lS). 

We require events to contain an T(1S) candidate by requiring that the two high momen- 
tum tracks in the event have an invariant mass in the range 9300 to 9600 MeV, consistent 
with the T(1S) mass. We make no cuts on the track quality variables for the lepton can- 
didate tracks, nor do we require them to satisfy lepton-identification criteria. The invari- 
ant mass requirement alone provides a nearly background free sample, and imposition of 
lepton-identification criteria only leads to larger systematic uncertainties and reduced signal 
efficiency without much improvement in signal quality. 

We require events to have three or four showers in the calorimeter, each of which has 
£' > 30 MeV, and is not matched to any charged track. Two of these showers must form 
an invariant mass within three standard deviations {3a) of the known 7r° mass. These 
candidates are kinematically constrained to the known 7r° mass, in order to improve the 
momentum resolution of the vr'' candidates. In addition to the two showers that correspond 
to the 7r°, events must contain an isolated photon candidate, between 50 and 250 MeV in 
energy, that does not form an invariant mass within 8 MeV {1.5a) of the 7r° mass with any 
other shower. Furthermore we require that the polar angle {6) of the third shower satisfy 
I cos^l < 0.804, the angular region in which CLEO's energy resolution is best. We allow 
events to contain up to one additional shower in the range |cos0| < 0.804. In addition, 
we allow for the possibility of one spurious charged track candidate in addition to the four 
"signal" tracks. Such spurious tracks may arise from failures in pattern recognition or from 



3 



1160903-007 



10800|- 



10500 



> 10200 

Q) 



9900 



9600- 



9300 



BB Threshhold 



rib(2S) T(2S) 

- iT 



Xb(2P) hb(2P) 







D(1D) 



^bPP ) hb(1P) 



+o- ..+ ■■- 



e"^e , ^i"^ [i 



0) 



0,1,2 



++ 



+- 



1,2,3- 



FIG. 1: The spectrum of bottomonium states below open bottom threshold. Dilepton decays of the 
J^'-' = 1 states are denoted by the solid lines, while dipion transitions are denoted by dashed 
lines. The radiative decay T(3S)— >7X6i,2(2-P) is indicated by the dot-dashed line, and the decay 
X6i,2(2-P)^wT(lS) by the dotted line. 



delta rays. Spurious showers may arise from synchrotron radiation from the beams or as 
a result of random noise in the calorimeter. If a given event yields more than one candidate 
due to the presence of an additional shower or track, we choose the candidate for which the 
sum of energies of all final state particles is nearest M(T(3S)). 

Because there is no phase space for a pair of kaons for decays in which a T(1S) is present, 
we assume that the low momentum charged tracks are pions. The invariant mass of the 
7r"'"7r~7r° combination, plotted in Fig. El exhibits a clear enhancement at the mass of u, 
M{u) = 0.783 GeV 0. 

To complete the reconstruction of the full decay chain, T(3S)— s>7c<jT(lS)-^77r+7r^7r°£+£~, 
we require that the /d.o.f. be less than two for a kinematic fit of tt^h'tt^ constrained to 
the uj mass, and subsequently that the mass recoiling against the kinematically fitted uj 
candidate and the photon lies within tf^ MeV of M(T(1S)) = 9.460 GeV 0- 

The simplest explanation for the observed events is the decay sequence 
T(3S)— >7Xfei,2(2P), with Xbi^2(2P)^co'T(lS). The lowest mass Xf,j(2P) state, Xbo(2-P), lies 
below threshold for decay to a;T(lS). In principle a transition through the rii,{?)S) state (see 
Fig. is possible, but this state has never been observed, and, furthermore, the energy 
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FIG. 2: The 7r"''7r~7r invariant mass, for data events subject to final analysis cuts with the exception 
of the cut on the /d.o.f. of the kinematic fit to uj. The overlaid histogram shows signal Monte 
Carlo events (normalized to the same total number) and indicates the good reproduction by the 
Monte Carlo of the shape and location of the u: peak. 

of the photon in the decay T(3S)^7r7h(3S') is expected to be below the range of observed 
energies in the data. 

Backgrounds from ordinary udsc quark pair production are extremely small because of 
the presence of the T(1S)— decay in the signal sample. The only significant source of 
background T(1S) expected is known cascades from T(3S). The final state of 77r"'"7r~7r°T(lS) 
may be reached through: 

T(3S)^7Xfc(2P), X6(2P)^7T(2S), T(2S)^7r+7r-T(lS) 
or T(3S)-^7rVT(2S), T(2S)^7r+7r-T(lS). 

In the first case, the final state of interest can be produced by the addition of a spurious 
shower in the calorimeter. In the second case, it may be reached by loss of one photon from 
one of the neutral pions due to acceptance or energy threshold. 

Backgrounds produced through either of these two processes are removed by excluding 
events in which the mass recoiling against the two charged pions in the T(3S) reference 
frame is consistent {i.e., between 9.78 and 9.81 GeV) with the hypothesis that the tt+tt" 
system is recoiling against T(2S) in the process T(3S)^X + T(2S)— >X + 7r+7r^T(lS). 

Two other decay sequences can yield the final state of 77r+7r^7r°T(lS): the decay 
T(3S)^7r+7r-T(2S), with the T(2S) decaying either to 7r07r°T(lS), or to '^Xb{lP) followed 
by X6(1-P)~^7'^(1S). In each of these cases, however, the charged pions are too soft to 
produce spurious uj candidates for the signal decay chain. 
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In order to evaluate residual background, we generated a GEANT j8| Monte Carlo sample 
for the channel T(3S)^7Xfei,2(2P),Xw,2(2P)^7T(2S), T(2S)^7r+7r-T(lS), corresponding 
to 21.5 ± 3.1 million T(3S) decays, or 4.53 ± 0.65 times our data set. The uncertainty on 
the equivalent number of T(3S) decays is due to the error on the branching ratios needed 
to convert our number of generated events to the equivalent number of T(3S) decays. This 
Monte Carlo sample produced one event that satisfied our selection criteria. We therefore 
expect 0.22 ± 0.03 events due to this source. We also generated a Monte Carlo sample of 
T(3S)^7r°7r°T(2S), T(2S)-^7r+7r-T(lS), corresponding to 540t^J° million T(3S) decays, or 
llStll times our data set. From this sample, a total of nine events passed our selection. 
In our data set, we thus expect 0.08 ± 0.01 events due to this source. To account for the 
residual background, we subtract the expected contribution of 0.30 events from the observed 
yield. We conservatively set a systematic error of ±0.30 events due to this subtraction. 

To evaluate the signal detection efficiency e, we generated 150,000 Monte 
Carlo events for each of Xfei(2-P) and Xb2{'^P), proceeding through the sequence 
T(3S)— ^>7Xbi,2(2-P)— *-7a;T(lS)— >-77r+7r~7r°£"'"£~, and uniform angular distributions for the 
T(3S)— >7Xbi,2(2-P) and T(1S)— decays. The masses for all particles in the decay chain 
were taken from Ref. 0- 

The analysis cuts described above are applied to these samples, and we obtain 
e(xfei(2P)) = (6.81±0.07)% and e(xb2(2P)) = (6.23±0.06)%, including all selection criteria, 
acceptance and trigger efficiencies. We apply an additional relative systematic error of ^3% 
to the efficiency in order to account for the possibility that the T(1S) retains the initial 
polarization of the T(3S). 

In order to illustrate the purity of the signal, in Fig. |2l we present a scatter plot of the 
mass recoiling against the Lwy system versus the dilepton invariant mass for all events subject 
to all the cuts discussed above, except those on the variables plotted. The final E.y spectrum 
is shown in Fig. |3J The observed yield has possible contributions from both decay sequences 
involving Xfei(2P) and Xfc2(2P) intermediate states. 

To obtain branching fractions for the X62(2P) and Xf,i(2P) transitions, we perform a 
maximum likelihood fit of the spectrum. The expected photon spectra for T(3S) transi- 
tions to X62(2P) and X6i(2P) were obtained from the signal Monte Carlo samples, and the 
observed spectrum was then fit to normalized Monte Carlo lineshapes with intensities 
(or yields) for XbiC^P) and X62(2P) as the free parameters. We obtain yields of 32.6^1? and 
20.1^5;! events, respectively. These yields have statistical significances of 10.2ct and 5.2a, 
respectively, obtained by comparing the likelihood of our final fitted yield to those of fits 
with zero signal events. The histogram resulting from the best fit is shown superimposed 
on the data in Fig. HJ 

The expected background contribution of 0.30 events is subtracted from the fitted yield 
by assuming that it scales as the ratio of the individual yields to the total yield. 

We thus obtain the following product branching ratios, using the detection efficiency and 
number of T(3S) decays discussed above: 



i3(T(3S)-^7X6i(2P)) X B{xbi{2P)^uTilS)) 
xl3{u^7T+7T-n^) X i3(T(lS)^£+r) 
= i0.82to\l ± 0.06) X 10'^ and 
i3(T(3S)^7X62(2P)) X Bixb2i2P)^ujTilS)) 
xi3(u;-^7r+7r-7r°) x S(T(lS)^£+r) 
= (0.55+[]:}^ ± 0.04) X 10~^ 



(1) 



(2) 
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FIG. 3: Plot of recoil mass versus dilepton mass, for data events subject to the final set of cuts, 
with the exception of the cuts on the two variables plotted. The number of events represented by 
each square range from one to four. 



in which the first uncertainty is statistical and the second is systematic. 

The statistical error of nearly 20% is dominant. Systematic error contributions to the 
above product of branching ratios are the following: 2% uncertainty in the number of T(3S), 
1% per charged track (a total of 4%) for track finding, 5% for 7^^ reconstruction, 2% for 
radiative 7 reconstruction, 1% for Monte Carlo statistics, lo% for the assumption of uniform 
T(1S)— angular distribution, and 1% for background subtraction. These contributions, 
added in quadrature, result in an overall relative systematic error of 

Using the present world average branching fractions for T(3S) ^7x61^2 (2-P), 
T(1S)— (taken to be twice that of T(1S)— and u—>-7t~^71~71^, we obtain 

B{x,i{2P)^coT{lS)) = (1.63l°:il°:J^)% and (3) 
i3te2(2P)-^T(lS)) = (1.10l°:il°:iJ)%, (4) 

The systematic errors include the additional uncertainty on the branching ratios for 
T(3S)^7Xw,2(2P), cj^7r+7r-7r° and T(lS)^£+£-, which contribute at the level of 5.9%. 

We may also calculate the ratio of Xfe2(2-P) to X6i(2P) branching ratios, for which several of 
the systematic errors discussed above cancel. We obtain this through a maximum likelihood 
fit to the spectrum, in which the two free parameters are the sum of yields and the ratio 
of Xb2(2P) to Xbii'^P) yields. When this fit is performed, we obtain a sum of yields equal 



7 



10 



>8 

c\i 

c/34 
c 

CD 

i£i2 



1160903-005 

I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 1 1 1 I I I I I I I I I I I I I I I I I 



I I I I h t IT^I 




j 

j . 

i : 

J 

.— j 
.j i 



Kinematically 
Forbidden 
Region 



-'1 r I I I I I I I I '-t--! — I- 



I ' ' ' ' I ' ' ' ' I 



I I I 





50 60 70 80 90 100 110 120 130 140 150 

E(y) IVIeV 



FIG. 4: Fitted photon energy spectrum for the final selection of events. The solid histogram shows 
contributions for both Xbii'^P) and Xb2(2-P), while the dotted and dashed histograms show the 
individual XbiC^P) and X62(2-P) contributions, repsectively. The dotted line indicates the region 
above which 7 energies are disallowed for lack of phase space. 



to 52.4lg;9 and a ratio of 0.62tQ;2o- In order to convert the ratio of yields to the ratio of 
branching ratios, we multiply the yield ratio by a factor of e(xw(2P))/e(xfe2(2P)). 

Bixb2i2P)^ujT{lS)))/Bixbii2P)^ujTilS)) (5) 

= 0.67l°i°. (6) 

The only systematic errors which do not cancel in this ratio are the small uncertainties in 
efficiency and T(3S)— >7X6i,2(2-P) branching ratios. These are entirely negligible compared 
to the output statistical error from the maximum likelihood fit. 

In Ref. j^, Voloshin predicts on the basis of S-wave phase space fac- 
tors, for El * El * El gluon configurations expected by the multipole expansion 
model 0, that r(x62(2P)^a;T(lS))/r(xw(2P)^cuT(lS)) ^ 1.4. The ratio of full 
widths r(x52(2P))/r(x5i(2P)) lies in the range of 1.25-1.5, using world average 
measurements of i3(xfei,2f 2P)— i>7T(lS, 2S)) and theoretical predictions for the rates 
r(Xfei,2(2P)— i>7T(lS, 2S)) jS]. Thus the branching ratios i3(xfei,2(2P)—>co'T(lS)) are expected 
to be approximately equal. Our measurement is in agreement with this expectation. 

In summary, we have made the observation of the first hadronic decay modes involving 
non-vector bb states using a sample of (5.81 ± 0.12) x 10^T(3S) decays collected by CLEO 
III. The discovery comes in the phase-space limited channels Xfei,2(2P)— >u;T(lS). We find 
that the ratio of the measured branching ratios for the two transitions are in agreement with 
theoretical expectations based on S-wave phase space factors for multipole expansions. 
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